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Nematic elastomers with a locked-in anisotropy direction exhibit semi-soft elastic response char- 
acterized by a plateau in the stress-strain curve in which stress does not change with strain. We 
calculate the global phase diagram for a minimal model, which is equivalent to one describing a 
nematic in crossed electric and magnetic fields, and show that semi-soft behavior is associated with 
a broken symmetry biaxial phase and that it persists well into the supercritical regime. We also 
consider generalizations beyond the minimal model and find similar results. 
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Nematic elastomers (NEs) [l[ are remarkable materials 
that combine the elastic properties of rubber with the ori- 
entational properties of nematic liquid crystals. An ideal 
uniaxial nematic elastomer is produced when an isotropic 
rubber, formed by crosslinking a polymer with nemato- 
genic mesogens, undergoes a transition to the nematic 
phase in which it spontaneously stretches along one di- 
rection (the z-direction) and contracts along the other 
two while its nematic mesogens align on average along 
the stretch direction. This ideal nematic phase exhibits 
soft-elasticity 0, 0] - a consequence of Goldstone modes 
arising from the breaking of the continuous rotational 
symmetry of the isotropic phase Soft elasticity is 
characterized by the vanishing of the elastic modulus C5 
measuring the energy associated with shears u^z and Uy^ 
in planes containing the anisotropy axis and by a stress- 
strain curve for strains Uxx (or Uyy) and stresses axx (or 
(Jyy) perpendicular to the anisotropy axis in which strains 
up to a critical value are produced at zero stress as shown 
in Fig. Ha). 

Monodomain samples cannot be produced without 
locking in a preferred anisotropy direction, usually by the 
Kiipfer-Finkelmann (KF) procedure Q in which a first 
crosslinking in the absence of uniaxial stress is followed 
by second one with stress. This process introduces a me- 
chanical aligning field h, analogous to an external electric 
or magnetic field, and lifts the value of the elastic mod- 
ulus C5 from zero. Thus, nematic elastomers prepared 
in this way are simply uniaxial solids with a linear stress 
strain relation at small strain. For fields h that are not 
too large, however, they are predicted to exhibit semi- 
soft elasticity [l|, 0] in which the nonlinear stress-strain 
curve exhibits a flat plateau at finite stress as shown in 
Fig. [Ija). Measured stress-strain curves in appropriately 
prepared samples unambiguously exhibit the character- 
istic semi-soft plateau 0, Q • 

The Goldstone argument for soft response predicts 
C5 = in the nematic phase, making reasonable conjec- 
tures that C5 should remain small at finite h when semi- 
soft response is expected and that semi-soft response 
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FIG. 1: (a): Soft (full line) and semi-soft (dashed and dotted 
lines) stress-strain curves at f = 0.08 with /i = 0, 0.8/ic, 4/ic, 
respectively, (b) Semi-soft curve of function of K^x 

at f = 0.08 and h = 2hc, where we have set v = w. 



might not exist at all in the supercritical regime [9| be- 
yond the mechanical critical point (with h = he) ter- 
minating the paranematic(PiV)-nematic(iV) coexistence 
line |10j| . There is now strong evidence [ll|, II4I that sam- 
ples prepared with the KF technique are supercritical. 
In addition, C5 measured in linearized rheological exper- 
iments is not particularly small These results have 
caused some to doubt the interpretation of the measured 
stress-strain plateau in terms of semi-soft response [13] ■ 

The purpose of this paper is to clarify the nature of 
semi-soft response. We consider the simplest or mini- 
mal model, which is formally equivalent to the Maier- 
Saupe-de-Gennes model Q for nematic liquid crystals, 
that exhibits this response. We derive the global mean- 
field phase diagram [Fig. [2] for this model. We show 
that semi-soft response is associated with biaxial phases 
that spontaneously break rotational symmetry, and we 
unambiguously establish that semi-soft response exists 
well into the supercritical regime. Figure [T] shows calcu- 
lated stress-strain curves for h = 0.8/ic and h = Ahc that 
clearly exhibit semi-soft behavior both for h < he and in 
the supercritical regime with h > he- Our minimal model 
provides a robust description of semi-soft response. We 
will, however, briefly discuss changes in this response that 
extensions of the minimal model can bring about. 

An elastomer is characterized by an equilibrium ref- 
erence configuration, which we refer to as a reference 
space Sb,, with mass points at positions x. Upon dis- 
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FIG. 2: Phase diagrams (a) in the h-a^x-f space showing the 
Sy and Sz {Sx hidden)CC and the Dx, Dy, and Dz DC 
surfaces along with the tricritical points tx, ty, tz and (b)in 
the f-h plane {d-^x = 0) showing the first-order uniaxial PN- 
N coexistence line NcCz, the mechanical critical point cz, and 
the Sz surface terminated by the line ft{h) with respective 
first- and second-order segments Nctz and tzE meeting at 
the tricritical point tz- 



tortion of the elastomer, points x are mapped to points 
R(x) = X -I- u(x) in a target space St, where u(x) is 
the displacement variable. Elastic distortions that vary 
slowly on scales set by the distance between crosslinks are 
described by the Cauchy deformation tensor A with com- 
ponents Ajj = dRi/dxj. The usual Lagrangian strain 
tensor is then u — (A-^ • A — 6)/2, where 6 is the unit 
matrix. The orientational properties of nematic meso- 
gens in the elastomer are measured by the Maier-Saupe 
nematic tensor Q. 

A complete theory for nematic elastomers should treat 
both A and Q and couplings between them. However, 
effective theories, obtained by integrating out Qij, that 
depend only on u provide a full description of the me- 
chanical properties of NEs 0, [1^1 • In such theories, 
strains measure distortions relative to an isotropic ref- 
erence state, and the elastic free-energy density f{u) 
consists of an isotropic part fiso{u) and an anisotropic 
part fani{u,h) arising from the imprinting process H|- 
Equilibrium in the presence of an external second Piola- 
Kirchhoff (PK) stress axx is determined by minimization 
over u of the Gibbs free energy density g{u, h, a^x, r) = 

fisoiu, r)+/ani(M, ^)+/cxt( ), where fcxt{u, 

—CTxxUxx- In equilibrium the second PK stress satisfies 
an — df/duij. We will return later to the engineering 



or first PK stress tensor a(j ~ df /dAij = AikCTkj- 

We now define our minimal model. First, we impose 
the constraint Tru = 0, enforcing incompressibility at 
small but not large u, rather than the full nonlinear in- 
compressibility constraint det A — [det{S + 2u)Y^^ = 1 
that more correctly describes NEs, whose bulk moduli 
are generally orders of magnitude larger than their shear 
moduli. Our theory thus depends only on the symmetric- 



traceless components of u: (jjij = mj — \5ijUkk, and 
/ext = —'^xx't'xx- Second , we use the simplest anisotropy 
energy: /ani = —huzz —h4>zz that favors stretching 
along the z axis. Thus, our theory is formally equivalent 
to that for a nematic liquid crystal in crossed electric 
and magnetic fields, E = Egz and H — Hex, in which 
4>i] ^ Qij, h ^ ^AeE^, and axx ^ \xaH'^, where Ae 
and Xa are, respectively, the anisotropic parts of the di- 
electric tensor and the magnetic susceptibility, and Gq, 
a = x,y, z, are unit vectors along direction a. Finally, we 
choose the Landau-de-Gennes form jl^] for fiso- 



/i,o(0,r) = irTr02 



(1) 



where we assume w > and where r — a(T — T*) with T 
the temperature and T* the temperature at the metasta- 
bility limit of the PN phase. In the isotropic phase with 
^ = 0, r = 2^, where is the T-dependent shear mod- 
ulus. We will often express quantities in reduced form: 



Uij = {v/w)uij, r = rv/w , h = hv 



C's = C^v/w , and similarly for other elastic moduli. 

We begin our analysis of the global phase diagram [l^ 
with the axx = plane, which we will refer to as the Z- 
plane because the anisotropy field h favors uniaxial order 
along the z-axis. The h> half of this plane exhibits the 
familiar nematic clearing point Nc at (rjv,/ijv) = (j^,0) 
and the PN-N coexistence line terminating at the me- 
chanical critical point {fc,hc) — (|i j^)- Throughout 
the h > half-plane, there is prolate uniaxial order with 
— S{ninj — ^Sij) with 5' > and the Frank director 
n along e^. In the TV phase at /i = and r < rjy, n can 
point anywhere on the unit sphere. Negative h induces 
oblate rather than prolate uniaxial order along and 
5 = — S" < at high temperature. When h < is turned 
on for r < rjv at which nematic order exists at /i = 0, n 
aligns in the two-dimensional xy-plane. This creates a bi- 
axial environment and biaxial rather than uniaxial order. 
Since n can point anywhere in the xy-plane, the biaxial 
state at h < exhibits a spontaneously broken symme- 
try. There must be a transition along a line r = rt{h) 
between the high-temperature oblate uniaxial state and 
the low-temperature biaxial state, which exists through- 
out the Sz surface shown in Fig. [51 This transition is 
first order at small \h\ because the PN-N transition is 
first order at /i = and second order at larger \h\, and 
there is a tricritical point [l3| tz sX {ft,ht) = ij^, "lili) 
separating the two behaviors as shown in Fig. ^h). A 
continuum of biaxial states coexist on Sz- We will re- 
fer to such surfaces as CC surfaces and ones on which a 
discrete set of states coexist as DC surfaces. 

The full phase diagram reflects the symmetries of g. 
The X- and ^-directions are equivalent in /iso, and the 
axx — and the h = planes are symmetry equivalent. 
These planes are also equivalent (apart from stretching) 
to the vertical plane with axx = h, but with positive and 
negative directions interchanged. To see this, we note 
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that (f>zz + 4)xx = -<t>vy and h(t)zz + crxx4ixx = -h4>yy 
when h = Gxx ■ Thus the phase structure of the Z-plane 
is replicated in the X-plane {h = 0) and the F-plane 
{(^xx = h) with respective preferred uniaxial order along 
e^; and Gy, critical points cx and cy, biaxial coexistence 
surfaces Sx and S'y, and tricritical points tx and ty- 

To fill in the 31? phase diagram, we consider perturba- 
tions away from the X-, and Z-planes. Turning on 
axx converts the PN-N coexistence line into a DC sur- 
face Dz, on which two discrete in general biaxial phases 
coexist. Turning on axx near the Sz surface favors align- 
ment of the biaxial order along ex when axx > and 
along By when axx < 0. Thus axx is an ordering field for 
biaxial order whereas a linear combination of h and axx 
acts as a nonordering field. The topology of the phase 
diagram near tz is that of the Blume-Emery-Griffiths 
model [3] with DC surfaces Dx and Dy emerging from 
the first-order line Nctz terminating Sz- The Dx and 
Dy surfaces terminate, respectively, on the critical lines 
Nctx and Ncty in the X- and F-planes. The surfaces 
Dx, Dy, and Dz form a cone with vertex at Nc- 

Before considering the axx-Uxx stress-strain curve, it is 
useful to look more closely at elastic response in the vicin- 
ity of the .Z-plane and the nature of order in the F-plane. 
Throughout the h > Q Z-plane, the equilibrium state is 
prolate uniaxial with order parameter S = Sq, and thus 
strains = ^Sq = — 2m°^ = —2Uyy. We are primarily 
interested in shears in the xz-plane and the response to 
an imposed axx with no additional stress along z. In this 
case Suzz — Uzz — u^zz '^il^ relax to an imposed Suxx, and 
the free energy of harmonic deviations from equilibrium 
can be written as Sf — ^C3{6uxx)^ + ^C5{6uxz)^ ■ The 
modulus C3 gives the slope of axx versus 5uxx, and C5 
is measured in linearized rheology experiments 12, 19| . 
C3 and C5 are easily calculated as a function of r and h. 
In reduced units, the ordered pair (CsjCs) takes on the 
value (|, ^) just above Nc (f = f^), (|, 0) just below Nc 
(f = f^), (0, ji^) at the critical point, and j^) in the 
supercritical regime at {f,h) — {fc,2hc). We will mea- 
sure elastic distortions using Suij rather than the strain 
relative to the reference space S'j^ defined by the equi- 
librium configuration at any given T j20] . 

On the h > 0, F-plane, there is oblate uniaxial order 
aligned along the y-direction at high T and biaxial order 
at low T. A convenient representation of the tensor order 
parameter is 
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where S' > 0. The vector ff = (771,772) = ri{cos29,sm29) 
is the biaxial order parameter, which is nonzero on the 
Sy surface. We define the equilibrium values of S' 
and 77 in the biaxial phase to be Sq and 770, respec- 
tively. Energy in this phase is independent of the ro- 
tation angle 9. Away from the K-plane, /ani + /oxt = 



-i(/i + axx)S' + {axx - h)r]i. Thus, axx < h favors 
7yi > and axx > h favors 771 < 0, implying that 
fj = (7^0,0) (or 61 = 0) at axx = and 77 = (-770,0) 
(or 9 = at axx = h'^ ■ These considerations im- 
ply that the modulus C5 is zero at axx — because 
C, = d^f/dulz\u.^^o = {2^o)-'d^f/d9^\e^o = 0. 

We can now construct the axx-Uxx stress-strain curves. 
At axx ~ 0, Uxx — Uxx', as axx is increased from zero, 
Suxx grows with initial slope I/C3 until axx = at 



is" - ?/" 

3*^0 '^xx 

to a maximum of ^Sq 



rjo. At axx = h, further 



770 pro- 



which point, 5uxx 
increase of duxx 

duces a zero-energy rotation of ff to yield duxx = J'^'q 
u'^.j. — 770 cos 29 and a nonzero shear Uxz = V2 = Vo sin 29. 
The growth of 772 from zero is induced by the vanish- 
ing of C5 at axx ~ and its becoming negative for 
|'7i| < '70- Thus, the characteristic semi-soft plateau is 
a consequence of Cs's vanishing at axx — h and not at 
(Jxx — 0. Measurements of C5 at axx = do not pro- 
vide information about what happens at axx — h. For 
f^xx > h, duxx again grows with axx- Figure [T] shows 
stress-strain curves for different values of h. Thus, semi- 
soft response is associated with the Sy surface, which 
exists at r and h well into the supercritical regime. 

A Ward identity provides a rigorous basis for the above 
picture beyond mean- field theory, /iso('u) is invariant 
under rotations of u, i.e., under u —^ UuU~^ where U is 
any rotation matrix. Thus if /ani = — Tr/iit, where hij — 
hez^ezJ, fiUuU-^) = fiM - TihUuU-\ for any U, 
including one describing an infinitesimal rotation by 7 



about the y-axis with components Ui^ 



^£yij-f, where 



eijk is the Levi-Civita anti-symmetric tensor. Equating 
the term linear in 7 in f{UuU~^) to that of TrhUuU~^ 
yields the Ward identity 



(Jxz [Uz 



c) = {(^zz + h~ axx)ux 



(3) 



where aij = df /duij. This identity applies for any /iso, 
including ones with no compressibility constraint, so long 
as /ani is linear in u. In the semi-soft geometry axz — 
azz — but axx > 0. Thus, either Uxz = or axx — 
h for any nonzero Uxz- Equation ((3]) also gives C5 = 
crxz/uxz\u^,^o = [h- (yxx)/{uzz - Uxx) ^ \h - axx\/2rio 
implying that C5 — > as axx — > as long as 770 ^ 0. 

We have focussed on the effects of an external second 
PK stress axx- In physical experiments, the first PK (en- 
gineering) stress, fjfj = df /dh-ij = AikC^kj, or the Cauchy 
stress, a^j — cr/^A^^/ det A (as in [H, Q), is externally 
controlled. The a^^-Axx stress-strain curve is easily ob- 
tained from the axx-Uxx curve using ct^^ = AxxO'xx and 
Axx — Vl + '^Uxx- These two curves are similar, but the 
flat plateau in the a^x-Axx curve rises linearly with Axx 
as shown in Fig. [lib), and there is a unique value of Axx 
for each value of cr^^. Thus, the Sy surface in the r-h-axx 
phase diagram would open into a finite volume biaxial re- 
gion in the r-h-a^^ phase diagram with a particular value 
of 7/ at each point in it. The phase diagram in the h-a^^ 
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FIG. 3: Schematic phase diagrams in the h{or /i')-stress plane. 
The points Tp, C and Cp are, respectively, triple points, 
liquid-gas-like critical points, and critical endpoints. (a) Di- 
agrams for the minimal model, where all transitions are first 
order; (b) and (c) Phase diagrams for more general /ani or 
/ext in which the first-order line from Sy is replaced by a 
surface terminated by two second-order (dashed) lines or one 
first-order line and one second-order line, h' and a'^^ are, 
respectively, the generalized aligning field and generalized ex- 
ternal stress resulting from the more general /ani or /ext. 



plane for < r < rt \s similar to that in Fig. [3jb). 

We have ignored boundary conditions and random 
stress, both of which can modify stress-strain curves. 
When Frank elastic energies are ignored, detailed cal- 
culations of domain structure induced by boundary con- 
ditions reproduce soft and semi-soft response fiT] . Small 
isotropic randomness appears not to affect soft response, 
but large randomness does p^ . Our approach should 
serve as a basis for further study of randomness. 

We can now consider modifications of the minimal 
model. A simple modification is to replace the constraint 
Trti = with the real volume constraint det A = 1. This 
replacement does not change the validity of the Ward 
identity and the resulting phase diagram has the same 
structure as that for Trit = but with different bound- 
aries for the CC and DC surfaces. In particular, the 
mechanical critical point is at ("re, he) = (0.1279,0.0052) 
and the tricritical point is at {ft,ht) = (0.1900,0.0247). 
Other modifications of the minimal model replace /ani 
with nonlinear functions of Uzz- Modifications of this 
kind can spread the CC surface Sy into a finite volume 
or convert it to a Z?C surface, as shown in Fig. [31 If 
/ani = ~hu'j,^, two states coexist, whereas with other 
forms such as might arise in a hexagonal lattice, three 
or more discrete states might coexist. When Sy is a 
DC surface, rather than exhibiting a homogeneous rota- 
tion of the biaxial order parameter (if boundary condi- 
tions are ignored) in response to an imposed u^x^ samples 
will break up into disrete domains of the allowed states. 
In other words, their response to external stress will be 
martensitic [2^ rather than semi-soft. 

The neo-classical model [24] can also be discussed in 
our language. The free energy of this model is a function 
of A and Q. It consists of an isotropic part, invariant 
under simultaneous rotations of A and Q in the target 
space and under rotations of A in the reference space, 
and a semi-soft anisotropic energy [6| , which is effectively 



nonlinear in the strain, that breaks rotational symmetry 
in the reference space. The phase diagram of this model 
is similar to that of the minimal model in the space of 
r-h-a^^. In it, semi-soft behavior also persists above the 
mechanical critical point ^25 ]. 

In summary, we determined the complete phase dia- 
gram of nematic elastomers subject to an internal align- 
ing field and a perpendicular external stress. Our results 
underscore the validity of semi-softness in the interpre- 
tation of their remarkable stress-strain curves. 

This work was supported by NSF grant DMR 0404570 
and the NSF MRSEC under DMR 05-20020. 
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